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Abstract. Previous work has shown that the cooling of SN 1987A excludes a Dirac- 
neutrino mass greater than (3(20 keV) for u e , i/ M , or v T . We re-examine the emission 
of wrong-helicitv, Dirac neutrinos from SN 1987A, and conclude that due to neutrino 
degeneracy and additional emission processes ( N + N N + N + uu, w~ + p — ► n + vv) 
the effect of a Dirac neutrino on the cooling of SN 1987A has been underestimated. We 
believe that the limit that follows from the cooling of SN 1987A is better — probably much 
better— than lOkeV. This result is significant in light of the recent evidence for a 17keV 
mass eigenstate that mixes with the electron neutrino. 
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I. Introduction 


Very-weakly interacting particles can be produced in the core of a newly born, hot 
neutron star, carry away energy, and accelerate the initial cooling process. Their observable 
effect is the shortening of the duration of the neutrino burst associated with the early 
cooling phase — and of course, the neutrino burst from SN 1987A was detected by the 
Kamiokande II (KII) and Irivine-Michigan-Brookhaven (IMB) water Cherenkov detectors. 
The potential shortening of the neutrino burst associated with SN 1987A has been used 
to severely constrain the properties of axions, 1 “righthanded” neutrinos, and other weakly 
interacting particles. Our interest here is in Dirac neutrinos; several authors have argued 
that a Dirac mass for any of the three neutrinos in the range of 0(20 keV) to 0(300 keV) 
is excluded by SN 1987A. 2 

The key to the argument involves the additional helicity states a Dirac neutrino has. 
While a majorana neutrino has but two helicity states: negative helicity (i/_), which 
corresponds to the neutrino state, and positive helicity (i/+), which corresponds to the an- 
tineutrino state, a Dirac neutrino has four: two helicity states associated with the neutrino 
(i/_, i/ + ) and two associated with the antineutrino (P_, P+). For a massless neutrino the 
helicity states, the eigenstates of a freely propagating neutrino, coincide with the chiral- 
ity states (left and right), the weak-interaction eigenstates of the neutrino: Vi = V- and 
Vr — is+ (or £c|_). In this case the additional, “wrong-helicity” states of a Dirac neutrino 
have no interactions — and are irrelevant. 

The situation changes if the neutrino has mass: the chirality and helicity states no 
longer coincide. In the ultrarelativistic limit the projection of V- (7"+) onto vi (vr) is 
order unity, and these helicity states have ordinary weak interactions. On the other hand, 
the wrong-helicity states, v+ (f>_), have but a small projection, order m l/ /2E u , on to 
the chirality states vi (vr), and to a first approximation are sterile. Owing to their 
small projection onto the weak-interaction chirality states they can be produced through 
ordinary weak interactions (“spin-flip” production): u~ —* v+ or i/ + — » p_. Of course it is 
also possible that the wrong-helicity states have other, new interactions (e.g., righthanded 
interactions). We will not address that possibility here. 

(Once produced, wrong-helicity neutrinos do interact through their projections onto the 
proper-helicity states; however for m u < 300 keV, the mean free path for such interactions 
is large compared to the size of a neutron star. For m„ > 300 keV wrong-helicity neutrinos 
should become trapped like their proper-helicity counterparts, and for a sufficiently large 
mass, their effect on the cooling will be comparable to that of the proper-helicity state 
neutrinos. The mass at which trapping is sufficient to make a Dirac species “supernova 
safe” must be greater than 300 keV, and an accurate determination of this mass requires 
a careful treatment of wrong-helicity neutrino transport. This is a formidable task. For 
our purposes it suffices to say that the value of the “supernova safe 1 ’ mass must certainly 
be greater than 300 keV, the mass where trapping sets in, and that for m < 300 keV 
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wrong-helicity neutrinos once produced stream out.) 

The most detailed study of the effect of Dirac neutrinos on the cooling of SN 1987A is 
that of Gandhi and Burrows. 3 In numerical models of the early cooling of SN 1987A they 
included the cooling effect of wrong-helicity neutrinos produced by the spin-flip-scattering 
processes, v~+N — > u + + N and A r +P+ — > A r + i'_. For neutron-star models cooled by both 
proper- and wrong-helicity neutrinos they computed the flux of proper-helicity neutrinos 
and the response of the KII and IMB detectors to this flux. They concluded that the 
duration of the detected neutrino bursts exclude a Dirac mass greater than about 14keV. 
In fact, their mass limit was extremely conservative; the effect of 14keV Dirac neutrino 
was to reduce the burst duration expected to less than about 1 sec in either detector. Had 
one instead insisted that the neutrino burst duration expected be no shorter than half the 
duration of the actual burst, their limit would have been about 9keV. 

On the face of it their work seems to preclude a Dirac neutrino of mass 17keV for 
example. Of course 17keV is a very interesting mass since several /3-decay experiments 
have found evidence for a 17keV neutrino-mass eigenstates that mixes with the electron 
neutrino at the 1% level (sin 2 0 ~ 0.01). 4 Moreover, the absence of neutrinoless double-/? 
decay in several isotopes strongly suggests that the l7keV mass eigenstate is of the Dirac 
type. Unfortunately, Gandhi and Burrows 3 recently discovered a simple factor of four error 
in the rate they used for the emission of wrong-helicity neutrinos, which has the effect of 
doubling their mass limit-raising their original limit to 28keV (and the less conservative 
limit that one could derive from their results to about 18keV). The motivation for re- 
examining the emission of wrong-helicity neutrinos from SN 1987 A hardly needs to be 
mentioned. 

To summarize our results briefly, we find that due to a number of effects the volume 
emissivity ( erg cm -3 sec -1 ) of wrong-helicity neutrinos is at least as large as — and probably 
much larger than — that used originally by Gandhi and Burrows, implying that their orig- 
inal “conservative limit” of 14keV stands. In particular the production of wrong-helicity 
neutrinos due to nucleon-nucleon, neutrino-pair bremsstrahlung is at least as important as 
spin-flip-scattering production. Since the cores of neutron stars are on the verge of pion 
condensation negative pions are likely to be present in great numbers. If this is the case, 
the process n~ + p — ♦ n + vv is likely to be even more important than the bremsstrahlung 
process (although less certain since the pion density depends critically upon the equation 
of state). Finally and probably most importantly, if there is significant mixing between the 
massive Dirac-neutrino (greater than few 0.1%), then deep in core of the neutron star the 
mass Dirac-neutrino species is degenerate (with chemical potential /i„ > 200 MeV) rather 
than nondegenerate as previously assumed; when this fact is taken into account the rate 
for wrong-helicity neutrino emission increases by a factor of order (^/T) 4 ~ 10 4 . 

While it is premature to quote a definitive limit to the mass of a Dirac neutrino based 
upon the cooling of SN 1987A, it seems clear that when all of the additional effects discussed 
here are incorporated into detailed numerical models of the early phase of neutron-star 
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cooling the mass limit will be more stringent than 10 keV, probably more like 1 keV. 


II. Spin-flip-scattering Production of Wrong-helicity Neutrinos 

Nondegenerate neutrinos 

Positive-helicity neutrinos (and negative- helicity antineutrinos) are produced by the 
helicity-flip scattering processes i/_ -f N — » v+ + N and + N — ► + N, where N is a 

nucleon. The matrix-element squared for this process has been computed by Gaemers et 
al., 5 

I-MsfI 2 = 8G 2 F mlm 2 [(4 + 3c 2 a ) - (c\r - c 2 A )cos6 ] , (1) 

where | Ad g p | 2 has been summed over initial and final nucleon spins, 8 is the angle between 
the incoming and outgoing neutrinos, m is the nucleon mass, m „ is the Dirac-neutrino 
mass, Gp ~ 1.17 x 10 -5 GeV 2 is the Fermi constant, and cy{p) — (1 — 4sin 2 6w)j 2 ~ 0, 
c a(p) — 9a/2, cy(n ) = —1/2, c^(n) ~ —g A /2, and <7,4 = 1.26. Unless stated otherwise we 
work in units where h = kg = c = 1. 

The volume emissivity ( erg cm -3 sec -1 ) of wrong-helicity neutrinos is given by 


£SF = J 1 A^sf j 2 (27r) 4 ^(p 1 + ki -p 2 - 


d 3 p 2 


) 3 2£ 2 (2tt) 3 


x- 


d 3 k\ d 3 k 2 


r/i(l — f% ')fi/k 2 , 


(2) 


2ki(2wy 2k 2 (2v) 3 

where p t is the four momentum of the incoming ( i = l)/outgoing ( i = 2) nucleon, k, is the 
four momentum of the incoming (i = l)/outgoing (i = 2) neutrino, /„ = [exp(Jfci/T) + l] -1 
is the phase-space distribution function of the incoming neutrino, and /, = [exp (Ej/T - 
fj>i/T) + 1]” 1 are the phase-space distribution functions of the nucleons. Note that we have 
allowed for nucleon degeneracy, but we have assumed that neutrinos are nondegenerate 
(in the standard scenario, a good assumption for v ^ and z/ r , but not v t ). Shortly we will 
return to the important issue of neutrino degeneracy. 

Making three reasonable assumptions this 12-dimensional integral can be reduced to 
a single integral. They are: (1) nonrelativistic nucleons; (2) negligible three momentum 
carried by the neutrinos (compared to the nucleons); and (3) incoming and outgoing neu- 
trinos have the same energy (elastic scattering). The volume emissivity for the process 
V- + N — » i/+ + N can then be written as 


£s.f — 


G 2 F ml(c 2 y + 3 c 2 a ) Pv { s/2 (mT) 3 / 2 d 
47 r j 7 r 2 dy 


j: 


y/udll 


■y + 1 


( 3 ) 


where y = (p i — m)/T and — 77r 2 T 4 /240 is the energy density in thermal neutrinos. In 
the nondegenerate limit, y <C — 1, the final term reduces to 71 Taking the nondegenerate 
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limit and taking into account the antineutrino process as well, the total volume emissivity 
becomes 

£ SF = 1.2 x 10 32 mi OO p 14 Tj O [0.9 + 0.2X„] erg cm -3 sec -1 , (4) 

where m \ oo is the neutrino mass in units of 100 keV, pu is the total mass density in units 
of 10 14 gem -3 , Tio is the temperature in units of lOMeV, and X n is the neutron fraction 
(the proton fraction X p = 1 - X n ). For arbitrary nucleon degeneracy, 

7G F m 2 m 3 / 2 T 11 / 2 . . . T t w t* \ 

£sf = 2 9 / 2 15 7r [l-2F(yj>) + 1.4J(y n )] , (5a) 

~ 2.6 x 10 31 m 2 00 (m/0.94 GeV) 3 ^ 2 T^ 2 (1.2/(y p ) + 1.4/(y„)j erg cm -3 sec -1 , (56) 

where I(y ) = ( d/dy ) f 0 °° y/udu/[exp(u — y) + 1] and we have displayed explicitly the de- 
pendence upon the nucleon mass m because the effective nucleon mass in nuclear matter is 
expected to be reduced by a factor of order 1/2. In the nondegenerate limit £sf does not 
depend upon the value of the nucleon mass, cf. Eq. (3); in the degenerate limit I(y ) oc y 1 / 2 
and y oc m~ 1 , so that £sf oc m 1 / 2 . The conditions at the core of the neutron star are ex- 
pected to be closer to nondegenerate than degenerate and £sf should be insensitive to the 
effective value of the nucleon mass. 

The following is a simple fit to I(y) that is accurate to better than 12% (typically 
accurate to a few %): 


ifn(y) -1 


2e - * 1 

y/* yfi + |y| 


1 

8(1 + |y|) 3/2 


Before going on to consider neutrino degeneracy, we should compare our expression for 
£sf with that used by Burrows and Gandhi. 3 After correcting the spin-flip cross section in 
their paper for the errant factor of four, we find that our volume emissivity (for X„ = 1/2) 
is a factor of about 1.9 larger than theirs. Most of the difference traces to one fact: We 
use |c j4 (n,p)| = <m/ 2 = 0.63 and they use 0.5. We believe that g A j 2 is the appropriate 
value to use. 6 The remaining discrepancy appears to involve round off (the spin-flip cross 
section used by Gandhi and Burrows is only given to one significant figure). 

Neutrino degeneracy 

At the core of a newly born hot neutron star, where the density is several times 
nuclear density and the temperature is of order 30MeV, electrons are highly degenerate 
with p e ~ Pr(e) ~ 240(A p pi 4 ) 1 / 3 MeV. On the other hand, neutrons and protons are only 
semi-degenerate, with p n , p p ~ 0(30 MeV). From this it follows that electron neutrinos are 
highly degenerate, as /5-equilibrium ( n + v e <-> p + e~) enforces: p„ e = p e + p p — p n ~ p e . In 
the absence of interactions that interconvert neutrinos of different flavor, p and r neutrinos 
should be nondegenerate. Since we know that electron-neutrino mass is less than about 10 
eV (more precisely, the mass of the dominant mass eigenstate associated with u e ), only the 
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production of wrong-helicity /j. and r neutrinos is of interest in setting a Dirac mass limit. 
In the absence of flavor-changing interactions they should be nondegenerate, justifying the 
previous assumption of nondegenerate neutrinos. 

However, mixing changes the story. If there is mixing between i/ e and v^, v T , then n 
and r neutrinos can become degenerate through u e <— > and i/ e <— ► v T oscillations. Since 
the matter-oscillation length of a neutrino is much less than its mean free path between 
weak interactions, the probability that an electron neutrino “next interacts as a /r (or r) 
neutrino” is simply sin 2 20 m /2, where 9 m is the mixing angle in matter between u e and 
(or v T ). The matter mixing angle is related to the vacuum mixing angle by: 7 


sin 29 m ~ sin 28 0 min [1, Aq/V], 


A 0 /V ~ 0.05 


/ 6m 2 \ ZlOOMeV 
V 10 8 eV 2 ) V ^ 





A 0 = 2E~ ^ — 3 pn eV , 

where E„ is the neutrino energy, 8m 2 is the difference of the mass squared between the 
two mass eigenstates, and V is the (weak-interaction) energy difference associated with the 
interaction of electron neutrinos, and /i and r neutrinos with the background (neutrons, 
protons, electrons, positrons, and neutrinos in the core of the neutron star). (Note, by 
using the min[l, Ao/V r ] we have not allowed for the possibility of resonant conversion; it 
would only further enhance neutrino mixing as sin 2 2 8 m — ► 1 in this case.) 

If neutrino mixing is effective, then chemical equilibrium will be established between 
the two (or three) neutrino species: //„ (and ji Vr ) = ^i l/c . One can estimate the rate at 
which neutrino mixing populates a degenerate sea of fi and r neutrinos: r sin 2 29 m /2 , 
where T v ~ n^G\E 2 v j n characterizes the rate at which neutrinos scatter off nucleons (the 
dominant neutrino interaction). Assuming Aq/V ^ 1, we estimate that the time required 
to populate a degenerate sea of fi (or r) neutrinos is 



That is, the mixing of a/ior r neutrino of mass of order lOkeV with the electron neutrino 
at the 1% level is sufficient to very rapidly populate a degenerate sea of /i or r neutrinos 
in the core of a hot neutron star. (Rapid here means much less than the cooling time of 
the neutron star: r < 1 sec). Note that neutrino oscillations mix flavors, but not helicity 
states, and so the degenerate sea of massive Dirac neutrinos filled by neutrino oscillations 
are proper-helicity neutrinos. The wrong-helicity neutrinos must still be produced by 
spin-flip processes. 

Degeneracy of /.i (and/or r) neutrinos will of course modify the “chemical composition” 
of the neutron star. While chemical equilibrium enforces fjL Ve = fi v = ii Vt and fj, Ve — 
Up ~ f*n + I'e, charge conservation, n c = 7T p , and lepton-number conservation, n t — + 

n Vt — n 0t + n Vii — n* + n Vr — riy r = n e — n € + 3(n l/e — n^> c ) — const, must also be observed. 
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Together these conditions determine all the chemical potentials. Qualitatively we can 
see that additional protons will have to decay to supply the additional neutrinos in the 
degenerate p and r seas: this will increase p n and decrease p p , p e , and p u . Since n oc p 
for a highly degenerate species, the neutrino and electron chemical potentials should not 
change by a large amount. However, a careful treatment of the effect of a massive Dirac 
neutrino that mixes with the electron neutrino on the initial cooling of a hot neutron star 
must take this into account. For our purposes we will assume that /i e ~ p v ~ 200 MeV. 

If the massive Dirac neutrinos are degenerate then our calculation of Csf must be 
revised: The neutrino distribution function must be changed to f v = [exp(E u /T — p„/T)-\- 
l] -1 . (No blocking factor need be added for the final-state neutrino since it is a wrong- 
helicity neutrino.) Making the same approximations as before, Eq. (3) is unchanged; 
however, the energy density of the neutrino p„ is now 


T 4 [°° u 3 du 

Pv ~2 ^ J 0 e^~y + 1 



(for p > T), 


where y — p v jT . In the highly degenerate limit, the energy density of neutrinos is much 
larger, order p 4 rather than order T 4 — there are more neutrinos and they have higher 
energies — and the volume emissivity is increased relative to the previous result, cf. Eq. 
(3), by a factor of 

10V/20OMeV)7T,V 

In the highly degenerate limit the process involving antineutrinos is severely suppressed 
as p , > = — /x„, implying that 


Pv 


t 4 r°° u 3 du 
2^2 J Q e “ +y + 1 


< Pv 


Bringing everything together, if we allow the neutrino sea associated with the massive 
Dirac neutrino to be degenerate, 


Glmlm 3 l 2 T u l 2 


£SF = 


2 5 /2 . 


7T° 


[1-2 1(y P ) + 1.4/(y n )] | J u 3 du 


+ 


1 


+ 1 e u +y + 1 


(6a) 


— ► 9.2 x 10 34 m 2 00 (m/0.94 GeV) 3/2 (/x^/200MeV) 4 T ? 0 /2 

x[1.2 1(y p ) + lAI(y n )] erg cm -3 sec -1 (for > T), (66) 

Note that the effect of neutrino degeneracy is always to increase £sf, because (p„ + pp) 
achieves its minimum for p u = 0. 

Provided that the massive Dirac neutrino is degenerate, the volume emissivity of wrong- 
helicity neutrinos is increased by a factor of 0( 10 4 ), which naively should improve the mass 
limit by a factor of 100. Holding the mixing angle fixed at 1%, we see that the timescale 
for populating the degenerate sea becomes of order 10 _1 sec for a mass of 1 keV, in which 
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case there is barely enough time to populate the degenerate sea of massive Dirac neutrinos. 
The SN 1987A mass limit clearly depends upon the mixing angle, and for 1% mixing it 
should be about 1 keV. 


III. Bremsstrahlung-pair Production of Wrong-helicity Neutrinos 

Wrong- heli city neutrinos can be produced through another spin-flip process: nucleon- 
nucleon, neutrino-pair bremsstrahlung, N N — » N + N + i/p, where vv — v-v _ or v+i>+ 
and N is a nucleon. The rate for this process can be found by using the matrix-element 
squared calculated by Friman and Maxwell 8 and the phase-space volume calculated by 
Brinkmann and Turner. 9 

The volume emissivity for N + N^>N + N + v+ + v+ is given by 


£BREM = 


J S\M\ 2 (2tt) 4 6 4 


(Pi + P2 ~ P3 ~ Pi — <]l ~ < 72 ) 


fp 1 

X 2E x (2ny 


d 3 p 4 d 3 q x d 3 q 2 
2 j F 4 (2tt) 3 2oj 1 (27t) 3 2cu 2 (27t) 3 


/l/ 2 (l — /3)(1 ~ / 4 V 1 , 


(7) 


where p x are the four momenta of the nucleons, /, = [exp (E{/T — pi/T) + l] -1 are 
the nucleon phase-space distribution functions, q x 2 are the four momenta of the neu- 
trino/antineutrino, u )\ t 2 are the energies of the neutrino/antineutrino, ui\ is the energy of 
the wrong-helicity neutrino, S is the symmetry factor (a factor of 1/2! for any pair of 
identical particles in the initial or final state), and the matrix- element squared is to be 
summed over initial and final nucleon spins. To begin we will assume that neutrinos are 
nondegenerate. 

Friman and Maxwell 8 have calculated the matrix element for the ordinary (no spin flip) 
pair-production bremsstrahlung process in the one-pion exchange approximation; their 
matrix element can be used to obtain the matrix element for the spin-flip process by 
multiplying by m u /2uji. Doing so and pulling out the only factor in the matrix element 
that depends upon the neutrino energies the desired matrix element can be written as 


\M\ 2 = I-MfmI 2 



( 8 ) 


where u = + ^2 is total energy carried off by the neutrino and antineutrino. 

There are actually three different bremsstrahlung processes: neutron-neutron, proton- 
proton, and neutron-proton. The matrix elements for the first two are the same. Calculat- 
ing the matrix-element squared is a tedious process involving the square of the sum of eight 
different diagrams (four direct and four exchange). Friman and Maxwell 8 have computed 
the square of the sum of the direct diagrams and the square of the sum of the exchange 
diagrams, but not the interference term. From previous experience with nucleon- nucleon, 
axion bremsstrahlung, 9 for which the matrix clement has a similar structure, we know that 
in the nondegenerate limit the interference term is very small, while in the degenerate limit 
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the contribution of the interference term increases |-M| 2 by about 50% (over the incoher- 
ent sum of the direct and exchange terms). Based on this we will ignore the interference 
terms (thereby likely underestimating the matrix element in the degenerate limit). The 
two matrix-elements squared (and summed over nucleon spins) axe 

I^FM(nn, pp)| 2 = 2 • 2 10 G 2 F g 2 A f i l-^FM( n P)| 2 = 3 ' 2™G 2 F g\f* » 

where G F = 1.17 x 10~ 5 GeV~ 2 is the Fermi constant, g A ~ 1.26 is the axial vector 
coupling constant, / ~ 1.1 is the pion-nucleon coupling, m is the nucleon mass, and 
m n ~ 0.135 GeV is the pion mass. Note: (1) we have already factored out the neutrino- 
energy dependence from the matrix-element squared; and (2) as presented, the Friman- 
Maxwell matrix-element squared must be multiplied by a factor of (2m) 4 because of their 
nucleon-spinor normalization convention. 

Now the phase-space volume integration. Brinkmann and Turner 9 have evaluated the 
five-particle phase-space volume element for nucleon-nucleon, axion bremsstrahlung for 
arbitrary nucleon degeneracy with the following assumptions: (1) nonrelativistic nucleons; 
(2) negligible axion three-momentum (compared to that of the nucleons); and (3) constant 
matrix element. As we will see very shortly the six-particle phase-space integral needed here 
can be reduced to the very same phase-space integral. For the neutrino bremsstrahlung case 
the analogous assumptions are: (1) nonrelativistic nucleons; and (2) negligible neutrino 
three momenta. In the axion case the integral over the axion’s momentum in the expression 
for i a that is analogous to Eq. (6) can be reduced to a single integral over the axion’s 
energy: 

In the present circumstance the integral over the momenta of the neutrino pair can be 
reduced to an integral over the sum of their energies: 


f d 3 9i 

d 3 q2 cj] Cu’2 


2u.' 1 (2tt) 

3 2u;2(27r) 3 oj 2 

\2u>i / V47T 2 y 48 J 


u r 


2 duj. 


Thus, by simply multiplying the results of Brinkmann and Turner 9 for the phase-space 
integrals over pi, ■ • • ,p 4 and E a by a factor of m 2 /1927r 2 we can obtain the phase-space 
integration needed here. 

The axion phase-space integral can be expressed as 

' ' ' 2E t Q*)> ^E a (2^6 4 ( P , +P2 - P3 - Pi - «.)/,/,(! - M( 1 - M 


I 


m l l 2 T™l 2 expfvi + V2)/1407 t 13 / 2 (nondegenerate limit), 
m i/2ji3/2 /( yi ,y 2 ) (in general), 


( 9 ) 
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where yi = (m — m)/T, and I(yi,y 2 ) is a (different) dimensionless function that must be 
evaluated numerically (see below). In the nondegenerate limit, 


n, = 2 





where Uj is the number density of species i (i = neutron or proton). 

To begin, consider the nondegenerate limit, a reasonable approximation to the condi- 
tions that pertain. 9 In this circumstance the volume emissivity can be expressed as, 


£BREM = 


16 g ^/ 4 

1057T 1 */ 2 mi 


f(Xn)n 2 N , 


(10a) 


^BREM - 1-5 x 10 31 f(X n erg cm 3 sec x , (106) 

where n n is the total nucleon density, A' R is again the neutron fraction, the function 
f(X n ) = 0.5 + 2X n (l — X n ) varies between 0.5 (for X n — 0, 1) and 1,0 (A" n = 1/2), and 
we have included a factor of two to account for both the process where the neutrino has 
the wrong helicity and the one where the antineutrino has the wrong helicity. We can 
compare this energy-loss rate to that from neutrino-nucleon spin-flip scattering (taking 
the nondegenerate limit for both): 


^BREM ^ q 14 P 14 f(X n ) 

i S F " ' T// 2 0.9 + 0.2X„‘ 

At the core of the newly born neutron star where most of the emission of wrong-helicity 
neutrinos occurs pu ~ 4 — 10 and T\o ~ 3 — 10, and thus the bremsstrahlung process 
should be of comparable importance. 

In the general case the volume emissivity is 



£brem = 


160/ 4 y 4 CF^m 9 / 2 T 13 / 2 
157r 2 m\ 


[0.5 {/(yi,^) + /(y 2 ,y 2 )} + 3/(yi,y 2 )] , 


(12a) 


~ 2.4 x 10 35 (m/0.94 GeV) 1/2 m 2 00 T™ /2 

x [0.5 {/(yi,yi) + 7(y 2 ,y2)} + 3/(yi,y 2 )] erg cm -3 sec -1 , (126) 

where we have displayed explicitly the kinematical dependence upon the nucleon mass 
(i.e., we have not pulled out the m 4 factor associated with the pion-nucleon coupling, 
m 4 /m 4 n ). On the basis of the nonlinear-sigma model it has been argued that the ratio 
of the nucleon mass to the pion mass should not change significantly with density. 10 In 
the nondegenerate limit, £brem varies as m ~ 5//2 and would increase by a significant factor 
if the effective nucleon mass is half its vacuum value. In the degenerate limit £brem is 
independent of the effective nucleon mass. 
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Since the nucleons in a newly born, hot neutron star are closer to being nondegenerate 
than degenerate, £brem oc m -5 ^ 2 will increase by a factor of 0 ( 6 ) if the effective nucleon 
mass is half its vacuum value, while £sf oc m° does not change. Thus, if the effective 
nucleon mass is substantially smaller than its vacuum value, the numerical factor in Eq. 
(11) is closer to unity, implying that the bremsstrahlung process dominates the spin- flip 
scattering process. 

Brinkmann and Turner 9 give a simple fit to /(yi,J/ 2 ) that is accurate to better than 
25% for all values of y i and 3 / 2 : 

htiyuVi)' 1 = 2.39 x 10 5 + 0.25e~*' + 0.25e-» 2 ) 

+1.73 x 10 4 (1 + |y |) — 1/2 + 6.92 x 10 4 (1 + |y|)~ 3/2 + 1.73 x 10 4 (1 + |y|)~ 5/2 , (13) 

where y = (yi + y 2 )/ 2 . 

If one is interested in producing helicity-flipped electron neutrinos, the URCA process 
can be very important (n + p — * n + n + e + + and p + p — > n + p + e + + u+). [Note the 
process where an electron rather than a positron is produced in the final state is highly 
suppressed because of electron degeneracy: p e ~ 300 MeV.] The matrix element for this 
process is four times larger than that for the neutron-neutron or proton-proton process. 
However, we are interested in the production of wrong-handed // and r neutrinos. 

Finally, Grifols and Masso 11 have also calculated the volume emissivity due to the 
bremsstrahlung process in the nondegenerate limit. Our results in this limit are larger 
than theirs by about a factor of five. The difference traces to a number of factors. First, 
they forgot to take in account the exchange diagrams (about a factor of two); second, they 
forgot to account for both wrong-helicity neutrino and antineutrino emission; and third, 
they did not take into account the pp bremsstrahlung process, which is important since 
during the early cooling phases X n ~ X p ~ 1/2. 

Neutrino degeneracy 

In computing the rate for the bremsstrahlung process we have assumed that proper- 
helicity neutrinos are nondegenerate. In light of our discussion of the important effect 
of neutrino degeneracy upon the spin-flip-scattering process we should re-examine that 
assumption. 

The effect here is far less pronounced. If the neutrino seas are degenerate, then there 
will be a significant blocking factor that suppresses the emission a proper-helicity neutrino, 
but not a proper-helicity antineutrino. In the nondegenerate limit the two processes, 
N+N — > N+N+v + +v+ and N+N —> N+N + u^+P-, contribute equally to £brem> in the 
highly degenerate limit only the first of these will contribute, the second being suppressed 
by the degenerate sea of v_ ’s. The net effect is a reduction of the volume emissivity by a 
factor of two. However, in this limit the spin-flip-scattering process is enhanced so much 
that the bremsstrahlung process becomes subdominant and unimportant. 
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Pion-nucleon , neutrino -pair production 

The conditions at the core of a neutron star are very close to those where a pion 
condensate should form. Because of this, the abundance of negative pions may well be 
comparable or greater than that of nucleons. 12 Needless to say, the abundance of pions in 
the core depends critically upon the equation of state. If the pion abundance is large, then 
the process 7r _ + p — > n v+ + v+ may also be important — in fact it may be dominant. 13 

The matrix element squared for this process is 


\M\ 2 = 


32G 2 F f 2 m 2 mlu>2 


m 2 u> i 


{^(1 - k • qi)(l - k • q 2 ) + 0.5(1 - 0.5k • q x - 0.5k • q 2 ) J , 


where we have assumed that the nucleons are nonrelativistic, that the pions are relativistic, 
and summed over initial and final nucleon spins. The three momenta of the pion, wrong- 
helicity neutrino, and proper-helicity antineutrino are k, q 1? and q 2 respectively, and the 
energies of these particles are uj (= uq +u^), aq , and u; 2 respectively. The volume emissivity 
for this process is given by 


EnN — 2 


/ 


d 3 pi d 3 p 2 d?k d 3 qi d 3 q 2 
2£'i(27r) 3 2£ 2 (27r) 3 2 w(2tt) 3 2w x (2tt) 3 2u; 2 (27r) 3 


x(2ir) 4 \M\ 2 6\p x +k-p 2 -q i -q 2 )f 1 f n u 1 , 

_ 5 (g 2 A + 0.5 )f 2 G 2 F ml(m/m n f 3 

— - I ApTl^Tl at , 

tt 3 rrr 

~ 2.1 x 10 33 mj 00 (m/0.94 GeV)~ 2 (n n /n^)X p pl 4 T 3 0 erg cm -3 sec -1 , 


(14a) 

( 146 ) 


where we have assumed that nucleons are nondegenerate, that the pion phase-space dis- 
tribution function f n — exp(a — k/T), and included a factor of 2 to account for both the 
process where the neutrino has the wrong helicity and the one where the antineutrino has 
the wrong helicity. Note that e n N oc m -2 , so that it is a factor of four larger if the effective 
nucleon mass is half its vacuum value. 

Now compare this production process to the bremsstrahlung process; taking the non- 
degenerate limit for f brem . cf. Eq. (10b), we find 


^ 7 rN 

£brem 


~ 140 



rp~ 1/2 
1 10 


(15) 


Thus, if the number density of negative pions is comparable to that of nucleons, pion- 
nucleon, neutrino-pair production is even more important than the bremsstrahlung process. 
In this case, unless the massive Dirac neutrino is degenerate, the pion-nucleon process will 
dominate. 


Axions 


In passing we note that it has been assumed that the dominant axion emission process 
for a hot young neutron is nucleon-nucleon, axion bremsstrahlung. 1 If there are lots of 
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negative pions in the core then pion-axion conversion 7r _ + p — ► n -f a, the analog of the 
process just discussed, can likewise dominate the bremsstrahlung process. In this case the 
axion volume emissivity is 


, _ 30 pgl N T> 

£a — 00 

(16a) 

7 r 


~ 4.9 x 10 49 X p ( — ) p1 4 Ti 0 erg cm -3 sec -1 , 

\n-N / 

(16! 0 


where g a w is a combination of axion-proton and axion-neutron couplings which is of order 
m/if a/ N ). The ratio of this process to the usual axion bremsstrahlung process is about 
50(n 7r /n j /v)T 1 "^ 1 ^ 2 , and if the abundance of negative pions is comparable to nucleons this 
process will be the dominant one. If this is the case, then the upper limit to the axion mass 
derived from SN 1987A improves by almost an order of magnitude: from about 10 -3 eV 
to almost 10“’ 4 eV. 1 * 14 

IV, Discussion 

We have computed the volume emissivity of wrong-helicity neutrinos due to the spin- 
flip scattering process off nucleons and nucleon-nucleon bremsstrahlung, both for arbitrary 
nucleon degeneracy and nondegenerate neutrinos. The two processes are found to be of 
comparable importance at the core of a newly born hot neutron star. Relative to the 
volume emissivity used by Gandhi and Burrows 3 (corrected for the errant factor of four) 
the total volume emissivity is larger by about a factor of four, which should restore their 
original, very conservative mass limit of 14keV. 

We have also calculated the production of wrong-helicity neutrinos (and antineutrinos) 
due to the process 7 + p — > n + uv, and find that if the number density of negative pions is 
comparable to that of nucleons (as could occur for a core on the verge of pion condensation), 
this process dominates both spin-flip scatterings and bremsstrahlung by a large factor. If 
this were the case the mass limit would improve to of order a few keV. Because the pion 
abundance is very sensitive to the equation of state, it is difficult to argue convincingly 
that such a bound is rigorous. 

Perhaps the most important effect is that of neutrino degeneracy. Electron neu- 
trinos are certainly degenerate at the core of a neutron star (with chemical potential 
fi e ~ 300 MeV); for masses in the keV range and mixing with the electron-neutrino of 
order 1% the massive Dirac neutrino should also become degenerate. This has the effect of 
increasing the emission of wrong-helicity neutrinos due to the spin-flip-scattering process 
by a factor of order (/u/T) 4 ^ 10 4 relative to the nondegenerate rates previously used. 
Provided that the massive Dirac neutrino mixes sufficiently with the electron neutrino this 
should improve the mass limit to around 1 keV. Of course, the massive Dirac neutrino 
need not mix with the electron neutrino at all. 

A precise limit to a Dirac-neutrino mass based upon SN 1987A awaits incorporation of 
the effects discussed here into detailed numerical cooling models, work which is currently 
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in progress. 15 However, it seems clear that the limit obtained will be more stringent than 
lOkeV, and if the massive Dirac neutrino mixes with the electron neutrino at the 1% level 
or more probably as stringent as 1 keV. 
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